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A  new  analytical  approach  is  proposed  for  evaluating  the  in-plane  permeability  of  gas  diffusion  layers 
(GDLs)  of  proton  exchange  membrane  fuel  cells.  In  this  approach,  the  microstructure  of  carbon  papers  is 
modeled  as  a  combination  of  equally-sized,  equally-spaced  fibers  parallel  and  perpendicular  to  the  flow 
direction.  The  permeability  of  the  carbon  paper  is  then  estimated  by  a  blend  of  the  permeability  of  the 
two  groups.  Several  blending  techniques  are  investigated  to  find  an  optimum  blend  through  comparisons 
with  experimental  data  for  GDLs.  The  proposed  model  captures  the  trends  of  experimental  data  over 
the  entire  range  of  GDL  porosity.  In  addition,  a  compact  relationship  is  reported  that  predicts  the  in¬ 
plane  permeability  of  GDL  as  a  function  of  porosity  and  the  fiber  diameter.  A  blending  technique  is  also 
successfully  adopted  to  report  a  closed-form  relationship  for  in-plane  permeability  of  three-directional 
fibrous  materials. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Proton  exchange  membrane  fuel  cells  (PEMFCs)  have  shown  the 
potential  to  be  commercialized  as  power  sources  in  automotive, 
electronic,  portable,  and  stationary  applications  [  1  ].  The  membrane 
electrode  assembly  (MEA),  the  heart  of  a  PEMFC,  is  composed  of 
a  membrane  loaded  by  catalyst  layers  on  each  side  and  is  sand¬ 
wiched  between  two  porous  layers  named  gas  diffusion  layers 
(GDLs)  [1].  In  addition  to  mechanical  support  of  the  membrane, 
GDL  allows  transport  of  reactants,  products,  and  electrons  from  the 
bipolar  plate  towards  the  catalyst  layer  and  vice  versa.  Therefore, 
thermophysical  properties  of  GDLs  such  as  gas  and  water  perme¬ 
abilities,  thermal,  and  electrical  conductivities  affect  the  PEMFC 
performance  and  reliability.  Transport  of  reactants  in  the  through- 
plane  direction  is  a  diffusive  mechanism  rather  than  a  convective 
mechanism;  thus,  the  only  pressure  driven  mass  transfer  occurs 
in  the  in-plane  direction  [2].  Several  theoretical  and  experimental 
studies  have  shown  that  the  in-plane  permeability  of  GDLs  is  a  key 
parameter  in  optimization  of  PEMFCs  performance  [3-8]. 

The  applied  pressure  for  sealing  a  PEMFC  changes  the  GDL  thick¬ 
ness  and  its  porosity;  this  affects  the  transport  properties  of  the 
GDL  For  example,  the  gas  permeability  reduces  as  a  result  of  the 
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GDL  compression  while  the  thermal  and  electrical  conductivities 
increase.  As  such,  a  trade  off  exists  in  the  design  and  optimization 
process  of  MEAs  [9-11].  Therefore,  having  a  model  that  can  pre¬ 
dict  the  in-plane  permeability  of  GDLs  as  a  function  of  porosity  is 
valuable  in  the  design  of  high  performance  PEMFCs. 

There  are  few  experimental  studies  available  in  the  literature 
that  reported  the  permeability  of  GDLs  [12-16].  Ihonen  et  al.  [14] 
showed  that  a  reverse  relationship  exists  between  the  in-plane  per¬ 
meability  and  compression.  Recently,  Feser  et  al.  [2]  studied  the 
effects  of  compression  and  reported  the  in-plane  gas  permeability 
as  a  function  of  porosity  for  a  carbon  cloth,  a  non-woven  carbon 
fiber  GDL,  and  a  carbon  paper.  In  a  similar  work,  Gostick  et  al.  [16] 
measured  the  permeability  of  several  commercial  GDLs  under  var¬ 
ious  compressive  loads  and  reported  the  in-plane  permeability  as 
a  function  of  porosity.  More  details  of  the  experimental  studies  can 
be  found  in  Ref.  [17]. 

She  et  al.  [18]  employed  a  fractal  permeability  model  that 
accounted  for  the  microstructures  of  GDLs  in  terms  of  two  fractal 
dimensions  and  proposed  a  relationship  for  permeability  as  a  func¬ 
tion  of  tortuosity,  fractal  dimensions,  pore  area  fractal  dimensions, 
pore  size  distribution,  and  effective  porosity.  Their  model,  however, 
requires  several  geometrical  parameters  (e.g.,  fractal  parameters) 
that  should  be  known  beforehand.  Lattice  Boltzmann  simulations 
of  gas  flow  through  several  random  fibrous  structures  were  car¬ 
ried  out  by  Vandoormaal  and  Pharoah  [19]  over  the  porosity  range 
of  0.6  <£<0.8.  They  reported  numerical  results  for  different  fiber 
orientations  which  were  an  order  of  magnitude  different  in  a  con¬ 
stant  porosity.  Based  on  the  permeability  values  for  structures  with 
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Nomenclature 

d  fiber  diameter  (m) 

GDL  gas  diffusion  layer 

K  viscous  permeability  (m2) 

I<eq  equivalent  permeability  of  the  mixture  (m2 ) 

P  pressure  (Pa) 

S  distance  between  centers  of  adjacent  fibers  (m) 

UD  volume-averaged  superficial  velocity  (m  s-1 ) 

Greek  symbols 

e  porosity 

/x  fluid  viscosity  (N  s  nrr2 ) 

cp  solid  volume  fraction,  cp  =  1  -  e 

ip'  dimensionless  parameter  in  Eq.  (3),  cp'  =  7t/4(1  -  s) 

Subscript 

par  parallel  to  flow  direction 

norm  normal  to  flow  direction 

tot  total 

eq  equivalent 


mixed  orientations,  they  proposed  two  correlations  for  in-plane 
and  through-plane  permeabilities  of  the  considered  structures 
which  were  not  compared  with  experimental  data  of  actual  GDLs. 

Our  literature  review  indicates  the  need  for  a  general  model 
that  can  accurately  predict  the  in-plane  permeability  of  GDLs,  and 
its  trends  as  a  function  of  porosity  and  fibers  diameter.  Therefore, 
the  objectives  of  the  present  work  are  to: 

•  Develop  and  verify  a  comprehensive  analytical  model  that  can 
predict  the  in-plane  gas  permeability  of  GDLs  and  other  complex 
fibrous  structures  that  captures  the  trends  observed  in  experi¬ 
mental  data. 

•  Investigate  the  effects  of  relevant  geometrical  parameters 
involved  and  identify  the  controlling  parameters. 

In  this  study,  a  novel  approach  is  proposed  for  studying  gas  flow 
through  GDLs  where  the  porous  medium  is  assumed  as  a  com¬ 
bination  of  fibers  normal  and  parallel  to  the  flow  direction.  We 
then  propose  that  a  blend  of  the  normal  and  parallel  permeabili¬ 
ties  of  unidirectional  fibers  can  provide  an  estimate  for  the  in-plane 
permeability  of  GDLs.  The  fraction  of  fibers  in  each  direction  is  eval¬ 
uated  from  a  simplified  ordered  geometry  structure  representing 
the  microstructure  of  carbon  paper.  Several  blending  techniques 
are  compared  against  experimental  data  collected  from  different 
sources.  It  is  observed  that  the  volume  weighted  blend  is  a  better 


choice  for  the  in-plane  permeability  of  GDLs.  Moreover,  a  compact 
relationship  is  presented  for  determining  the  in-plane  permeability 
of  GDLs  as  a  function  of  porosity  and  fiber  diameter.  The  applica¬ 
bility  of  blending  technique  for  predicting  in-plane  permeability 
of  fibrous  materials  with  three-directional  microstructure  is  also 
investigated.  The  volume-weighted  resistivity  scheme  captures  the 
trend  of  experimental  data  for  permeability  of  hydrogels,  metallic 
rods,  metalfoams,  and  glass  wool. 


2.  Model  development 


Permeability,  K,  can  be  interpreted  as  a  measure  of  the  flow  con¬ 
ductance  of  the  porous  material  and  is  defined  using  Darcy  equation 
[20]: 


dP  ^ii 

-Tx  =  kUd 


(1) 


where  Ud  is  the  volume  averaged  velocity  through  porous  media,  pi 
is  viscosity  and  dP/dx  is  the  pressure  difference  between  two  points 
in  the  flow-field.  Darcy  equation  with  the  form  of  Eq.  (1)  is  valid 
for  incompressible,  steady,  constant  properties,  single-phase  (no¬ 
surface  tension  forces),  and  low  Reynolds  number  flows.  As  such, 
the  same  assumptions  are  considered  in  the  present  analysis. 

There  are  different  types  of  GDLs  including  carbon  paper  type, 
carbon  fiber  cloth,  and  wet-/dry-laid  papers  [19].  The  carbon  paper 
type  is  the  focus  of  the  present  study.  An  SEM  image  of  Toray  carbon 
paper  is  shown  in  Fig.  1 .  Gas  permeability  of  such  fibrous  structures 
depends  on  several  factors  including:  porosity,  fibers  size  and  dis¬ 
tribution.  In  carbon  papers,  the  axes  of  fibers  are  located  on  parallel 
planes  with  random  distribution  and  orientation,  see  Fig.  1.  There¬ 
fore,  the  geometry  shown  in  Fig.  2  can  be  assumed  as  an  idealized 
representation  of  the  GDLs.  The  porosity  of  the  ordered  structure 
is  related  to  the  distance  between  the  centers  of  adjacent  fibers,  S, 


Fig.  1.  SEM  image  of  Toray  90  carbon  paper:  (a)  top  view;  (b)  side  view. 
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Fig.  3.  Summary  of  the  blending  technique  concept  for  2D  fibrous  structures  (GDLs). 


and  the  fibers  diameter  d : 


„  71  d 

e=1-4S 


Fig.  4.  Comparison  of  Eq.  (3)  with  experimental  data  for  normal  permeability  of 
unidirectional  fibers  [25]. 


Following  Jackson  and  James  [21],  the  carbon  paper  is  modeled 
as  a  mixture  of  fibers  parallel  and  normal  to  flow  direction  with 
solid  volume  fractions  of  (ppar  and  <pnorm,  respectively;  for  the  geom¬ 
etry  shown  in  Fig.  2,  (ppar  =  (Pnorm  =  (ptot/2 .  As  such,  it  is  expected  that 
the  equivalent  permeability,  Keq{(p ),  is  related  to  the  permeability 
of  each  component;  the  concept  of  blending  technique  is  shown  in 
Fig.  3.  It  should  be  noted  that  there  is  no  concrete  rules  for  esti¬ 
mating  the  mixtures  permeability.  Blending  techniques  have  been 
successfully  employed  to  estimate  the  permeability  of  fibrous  mix¬ 
tures  such  as  hydrogels  [21],  fibers  with  different  sizes  [22],  and 
fibers  with  different  charges  [23,24].  Flowever,  to  our  best  knowl¬ 
edge,  the  application  of  blending  techniques  to  planar  structures 
such  as  GDL  is  novel. 

Different  blending  rules,  originally  developed  for  different 
applications,  are  rewritten  for  a  mixture  of  fibers  with  different 
orientations  and  are  listed  in  Table  1.  One  difference  among  the 
blending  schemes  is  the  way  the  interactions  between  a  fiber  and 
its  neighbours  are  considered.  The  other  difference  is  the  way 
each  fiber  category  affects  the  overall  pressure  drop.  The  volume- 
weighted  permeability  is  based  on  the  assumption  that  the  two 
types  of  fibers  occupied  separate,  parallel  phases  in  a  slab  of  fibrous 
material  [24].  In  other  words,  the  volume-weighted  permeability 
model  considers  each  fiber  category  as  parallel  flow  resistors.  The 
weighted  resistivity  schemes  assume  that  the  two  fiber  categories 
act  similar  to  flow  resistors  in  series,  i.e.,  the  resistivity  based  meth¬ 
ods  consider  the  fibrous  material  as  combination  of  various  layers 
in  series.  The  geometric  mean  scheme  is  a  purely  mathematical 
blend. 

To  estimate  the  mixture  permeability  from  using  the  equations 
listed  in  Table  1,  one  needs  to  know  the  permeability  of  the  fibers 
in  both  normal  and  parallel  directions. 


Table  1 

A  summary  of  various  blending  methods  investigated  for  estimation  of  in-plane 
permeability  of  planar  structures  (GDLs). 


Blending  model 

Relationship 

Volume-weighted  resistivity 

1  _  Vpar  1  .  i Pnorm  1 

IQq(<p)  <P  Kpar{(p )  <P  Knorm{<p) 

Unweighted  resistivity 

1  -  1  +  1 

Keq(<P)  Kpar(<Ppar)  Knorm[<Pnorm ) 

Volume  weighted  permeability 

Keq(V>)=  ^Kpar(<p)+ 

Geometric  mean 

Keq(<p)  -  (Kpar((p)yp-’iv(Kmrm(<p)y™iv 

< Ppar  —  (Pnorm  — 

2.  \ .  Normal  permeability 

Recently,  Tamayol  and  Bahrami  [25-27]  have  analytically  stud¬ 
ied  the  normal  and  parallel  permeabilities  of  ordered  arrangements 
of  unidirectional  fibers.  In  Ref.  [25],  assuming  a  parabolic  veloc¬ 
ity  profile  within  square  arrangements  of  equally-sized  fibers  and 
integrating  the  continuity  and  momentum  equations,  a  closed  form 
analytical  solution  has  been  reported  for  the  pressure  drop  and 
permeability: 


(  12(yV-l) 

' 2-g(s ) 

[  <p'  \fv' 

2 

18 

!  18  +  12(^-1)  ( 

tan"{vh)  +  \ 

4 

d\ 

V Vu  -v? 

/ 

¥  =  TpJ—r,  ■  «(£)  = 1  -274£  -  °-274 

(3) 

Eq.  (3)  is  compared  with  experimental  data  for  normal  perme¬ 
ability  of  unidirectional  fibers  [28-32]  in  Fig.  4.  It  can  be  seen  that 
the  model  of  [25]  is  in  a  reasonable  agreement  with  experimental 
data  over  the  entire  range  of  porosity.  As  such,  Eq.  (3)  is  used  to 
calculate  the  permeability  of  fibers  normal  to  flow  direction  in  the 
blending  models. 


2.2.  Parallel  permeability 


Tamayol  and  Bahrami  [26]  studied  steady,  incompressible,  and 
fully-developed  flow  parallel  to  axes  of  mono-disperse  fibers  in 
several  ordered  arrays  of  cylinders.  They  started  from  general  solu¬ 
tion  of  the  Poisson’s  equation  and  using  point  matching  technique 
reported  closed-form  analytical  solutions  for  velocity  distribution 
in  the  considered  geometries  including:  square,  staggered,  and 
hexagonal  arrays  of  fiber.  Their  solution  for  parallel  permeability 
of  square  fiber  arrangements  is: 


d2 

16  <p 


-1.479  -  In  <p  +  2(p 


2 


0.01 86<p4 


(4) 


where  <p  is  the  fiber  volume  fraction,  (p  =  1  -  s.  In  Fig.  5,  this  model 
is  compared  with  experimental  and  numerical  data  found  in  the 
literature  [26,33-35].  As  shown  in  Fig.  5,  Eq.  (4)  is  also  in  good 
agreement  with  the  collected  data  over  the  entire  range  of  porosity. 
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Fig.  5.  Comparison  of  Eq.  (4)  with  experimental  data  for  parallel  permeability  of 
unidirectional  fibers  [26]. 


Fig.  7.  Comparison  of  the  proposed  blending  model  with  the  experimental  data  for 
planar  structures;  experimental  data  from  Feser  et  al.  [2]  and  Gostick  et  al.  [16]  are 
for  various  types  of  GDLs. 


Therefore,  Eq.  (4)  is  employed  to  estimate  the  parallel  permeability 
in  the  blending  rules. 

2.3.  Comparison  of  different  blending  rules  for  planar  structures 
(GDLs) 


This  relationship  is  only  a  function  of  porosity  and  fiber  diam¬ 
eter.  Comparison  of  Eq.  (5)  with  experimental  data  is  presented  in 
the  following  sections. 

3.  Comparison  with  experimental  data 


The  blending  models  listed  in  Table  1  combined  with  normal 
and  parallel  permeabilities  presented  by  Eqs.  (3)  and  (4)  are  used 
to  calculate  the  permeability  of  the  geometry  shown  in  Fig.  2.  In 
the  considered  geometry,  50%  of  fibers  are  parallel  and  50%  of  fibers 
are  normal  to  flow  direction,  i.e.,  (ppar  =  (pnorm  =  (Ptot/2 •  The  calculated 
results  from  different  blending  rules  are  plotted  in  Fig.  6.  It  can 
be  seen  that  the  blending  rules  fall  between  normal  and  parallel 
permeabilities  of  square  arrays  of  cylinders  and  at  high  porosities, 
s  >  0.7,  there  is  a  small  difference  between  different  models. 

Our  analysis  showed  that  the  volume  weighted  permeability 
scheme  is  in  better  agreement  with  experimental  data  in  lower 
porosities,  £  <  0.7.  This  blending  model  is  written  in  the  following 
easy-to-use  form: 


I<{£,  d)  =  exp 


f  -12.95  + 13. 9£  \  2 
VI  +1.576:  -2.22e2J  d 


(5) 


Fig.  6.  Comparison  of  different  blending  models  with  the  bounds  for 

f Ppar  =  (Pnorm  =  <P  tot  1 2. . 


The  blending  model  is  compared  with  experimental  results  of 
Gostick  et  al.  [16]  and  Feser  et  al.  [2]  for  a  variety  of  carbon  paper 
GDLs  in  Fig.  7.  In  addition,  to  increase  the  porosity  range  of  experi¬ 
mental  data,  the  results  reported  by  Shi  and  Lee  [36]  for  composite 
fabrication  application  are  included.  The  composite  reinforcement 
mats  have  a  similar  microstructure  to  GDLs  [37].  The  permeability 
values  are  non-dimensionalized  using  the  fiber  diameters  reported 
in  the  abovementioned  studies. 

Fig.  7  shows  that  the  volume  weighted  permeability  method 
predicts  the  trends  of  experimental  data  over  a  wide  range  of  poros¬ 
ity. 

As  such,  it  can  be  concluded  that  a  parallel  network  comprised  of 
the  flow  resistors  for  fibers  parallel  and  normal  to  flow  directions, 
captures  the  overall  flow  resistance  accurately. 

Most  of  the  experimental  data  fall  between  the  normal  and  par¬ 
allel  permeabilities  of  the  square  arrangement  of  fibers.  Therefore, 
normal  and  parallel  permeability  of  unidirectional  fibers  can  be 
viewed  as  upper  and  lower  bounds  for  the  in-plane  gas  perme¬ 
ability  of  fibrous  porous  media  such  as  GDLs,  i.e.,  the  permeability 
of  fibers  with  mixed  orientations  is  bounded  by  the  two  limiting 
cases  when  all  of  the  fibers  are  oriented  either  parallel  or  normal 
to  the  flow  direction.  This  is  in  line  with  the  previous  observations 
of  Tomadakis  and  Robertson  [37]  and  Tamayol  and  Bahrami  [25]. 

Our  analysis  suggests  that  in-plane  permeability  is  proportional 
to  the  fibers  diameter  squared;  this  is  in  line  with  other  theoretical 
models  listed  in  Table  2.  The  experimental  data  for  various  materi¬ 
als  with  different  diameters,  non-dimensionalized  using  d2  in  Fig.  7, 
follows  our  proposed  relationship,  which  verifies  existence  of  this 
proportionality. 

4.  Comparison  with  other  existing  models 

In  Fig.  8,  the  present  model,  Eq.  (5),  and  the  collected  exper¬ 
imental  data  are  compared  against  the  correlations  reported  by 
Tomadakis  and  Sotirchos  [40]  and  Van  Doormaal  and  Pharoah  [19]; 
see  Table  2  for  more  details.  The  model  of  Tomadakis  and  Sotirchos 
(TS)  [38]  was  based  on  the  analogy  between  electrical  and  flow 
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Reference 

Relationships 

Remarks 

[38] 

/C=  e  2(  (£-£p)(“+2)  ^d2,  ep  =  0.11,  a  —  0.521 

8(ln  e)2  y(-l-spf[(a+-l)s-ep]2  J  P 

•  Based  on  the  analogy  between  electrical  and  flow  conductions 

•  Accurate  only  for  s  <  0.8 

•  Developed  for  over  lapping  fibers 

[19] 

K  =  0.065  f^d2 

•  Based  on  curve  fit  of  numerical  results 

•  Accurate  only  for  0.6  <  e  <  0.8 

conductions.  This  model  was  originally  developed  for  permeability 
of  randomly  distributed  overlapping  fibers  in  composite  reinforce¬ 
ments  [37].  It  can  be  seen  that  both  Eq.  (5)  and  the  model  of  [38] 
predict  the  trends  of  experimental  data  over  the  low  to  medium 
range  of  porosity.  However,  TS  model  [38]  overpredicts  the  data 
in  high  porosities,  £<0.8,  while  Eq.  (5)  is  in  agreement  with  the 
experimental  data  over  the  entire  range  of  porosity. 


5.  Permeability  of  three-directional  fibrous  structures 


To  further  investigate  application  of  blending  techniques  for 
estimating  permeability  of  fibrous  media  with  complex  but  non- 
planar  microstructure,  a  similar  analysis  has  been  conducted  for 
three-directional  [37]  fibrous  structures  with  random  distribution 
and  orientation  of  fibers  in  the  space.  Following  Jackson  and  James 
[21  ]  the  complex  geometry  of  is  modeled  with  the  simple  cubic  (SC) 
arrangement  shown  in  Fig.  9.  In  the  SC  structure,  1/3  of  fibers  are 
parallel  and  2/3  of  fibers  are  normal  to  flow  direction  see  Fig.  10, 
i.e.,  (ppar  =  (ptotl 3,  (Pnorm  =  2(Ptotl3 .  The  volume-weighted  resistivity 
scheme  is  in  a  reasonable  agreement  with  experimental  data  and 
is  written  in  the  following  compact  form: 


I<(s,  d)  =  exp 


/  -43.25  +  46.6s  \ 

U  +  10.56s-  10.5s2  )  d 


(6) 


In  Fig.  1 1 ,  the  proposed  model  is  compared  with  the  experimen¬ 
tal  values  of  permeability  of  hydrogels,  metallic  rods,  metalfoams, 
and  glass  wool  [39-44].  Moreover,  a  selection  of  other  existing 
models  [21,38]  is  included  in  the  figure.  It  can  be  seen  that  the 
blending  technique  is  the  only  model  that  captures  the  trend  of 
experimental  data  over  the  entire  range  of  porosity. 


6.  Conclusions 


A  novel  analytical  model  was  developed  for  the  in-plane  per¬ 
meability  of  gas  diffusion  layers  as  a  function  of  porosity  and  fiber 
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Fig.  8.  Comparison  of  present  model  with  other  existing  correlations  for  the  in-plane 
permeability  of  fiber  mats. 


Fig.  9.  Proposed  periodic  geometry  for  modeling  3D  (non-planar)  fibrous  structures. 


diameter.  In  this  approach,  the  porous  medium  was  considered  as 
a  mixture  of  fibers  parallel  and  normal  to  flow  directions.  Then, 
the  permeability  of  the  mixture  was  modeled  as  a  blend  of  the 
permeabilities  of  its  components.  The  normal  and  parallel  fiber 
permeabilities  were  evaluated  from  our  previous  analytical  stud¬ 
ies  for  flow  through  square  arrays  of  fibers.  To  find  the  fraction 
of  each  fibers  category,  an  ordered  equally-sized,  equally  spaced 
planar  microstructure  was  assumed.  The  highlights  of  the  present 
study  can  be  summarized  as: 


Fig.  10.  The  blending  technique  concept  for  3D  fibrous  structures. 
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Fig.  11.  Comparison  of  present  model  and  other  existing  correlations  with  experi¬ 
mental  data  for  the  in-plane  permeability  of  three-directional  fibrous  structures. 

•  Normal  and  parallel  permeability  of  the  fibers  provides  bounds 
for  the  permeability  GDLs,  see  Fig.  7. 

•  GDLs  can  be  treated  as  a  combination  of  fibers  with  different 
orientations.  Therefore,  the  proposed  approach  may  be  used  for 
predicting  other  transport  properties. 

•  Permeability  of  GDLs  is  directly  proportional  to  its  porosity  and 
the  fibers  diameter  squared. 

Employing  the  blending  rules,  a  compact  relationship  was  pro¬ 
posed  for  gas  in-plane  permeability  of  GDLs  which  was  a  function 
of  the  porosity  and  fibers  diameter.  The  model  can  be  used  to 
guide  the  design  and  optimization  of  PEMFCs,  and  can  be  readily 
implemented  into  fuel  cell  models  that  require  specification  of  the 
in-plane  gas  permeability  of  the  GDL. 

The  blending  approach  with  the  volume-weighted  resistivity 
scheme  captures  the  trends  of  values  for  in-plane  permeability  of 
three-directional  fibrous  materials  such  as  hydrogels,  metallic  rods, 
metalfoams,  and  glass  wool. 
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